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ABSTRACT 

Motivated  by  the  potential  of  oceanic  mesoscale  eddies  to  drive  intrinsic  low-frequency  variability,  this 
paper  examines  geostrophic  turbulence  in  the  frequency-wavenumber  domain.  Frequency-wavenumber 
spectra,  spectral  fluxes,  and  spectral  transfers  are  computed  from  an  idealized  two-layer  quasigeostrophic 
(QG)  turbulence  model,  a  realistic  high-resolution  global  ocean  general  circulation  model,  and  gridded 
satellite  altimeter  products.  In  the  idealized  QG  model,  energy  in  low  wavenumbers,  arising  from  non¬ 
linear  interactions  via  the  well-known  inverse  cascade,  is  associated  with  energy  in  low  frequencies  and 
vice  versa,  although  not  in  a  simple  way.  The  range  of  frequencies  that  are  highly  energized  and  engaged  in 
nonlinear  transfer  is  much  greater  than  the  range  of  highly  energized  and  engaged  wavenumbers.  Low- 
frequency,  low-wavenumber  energy  is  maintained  primarily  by  nonlinearities  in  the  QG  model,  with 
forcing  and  friction  playing  important  but  secondary  roles.  In  the  high-resolution  ocean  model,  non- 
linearities  also  generally  drive  kinetic  energy  to  low  frequencies  as  well  as  to  low  wavenumbers.  Impli¬ 
cations  for  the  maintenance  of  low-frequency  oceanic  variability  are  discussed.  The  cascade  of  surface 
kinetic  energy  to  low  frequencies  that  predominates  in  idealized  and  realistic  models  is  seen  in  some 
regions  of  the  gridded  altimeter  product,  but  not  in  others.  Exercises  conducted  with  the  general  circu¬ 
lation  model  suggest  that  the  spatial  and  temporal  filtering  inherent  in  the  construction  of  gridded  satellite 
altimeter  maps  may  contribute  to  the  discrepancies  between  the  direction  of  the  frequency  cascade  in 
models  versus  gridded  altimeter  maps  seen  in  some  regions.  Of  course,  another  potential  reason  for  the 
discrepancy  is  missing  physics  in  the  models  utilized  here. 


*  Naval  Research  Laboratory  Contribution  Number  NRL/JA/ 
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1.  Introduction 

Geostrophic  turbulence  and  its  progenitor  two- 
dimensional  turbulence  serve  as  valuable  paradigms  for 
atmospheric  and  oceanic  flows.  This  paper  addresses  the 
relative  contributions  of  forcing  versus  intrinsic  non- 
linearities  in  the  maintenance  of  low-frequency  variability 
in  oceanic  geostrophic  turbulence.  Our  main  diagnostics 
are  spectral  fluxes  and  spectral  transfers — long  utilized  in 
wavenumber  space  analyses  of  geostrophic  turbulence — 
computed  here  in  frequency-wavenumber  space. 
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The  role  of  forcing  versus  intrinsic  nonlinear  pro¬ 
cesses  in  the  maintenance  of  low-frequency  variability  is 
an  important  topic  in  climate  dynamics.  The  climate 
system  exhibits  variability  over  a  vast  range  of  time  scales. 
Huybers  and  Curry  (2006)  show,  for  instance,  that  the 
frequency  spectrum  of  surface  temperature  from  monthly 
out  to  100000-yr  time  scales  displays  power  at  all  fre¬ 
quencies.  Some  low-frequency  variability  is  due  to 
forcing.  Hasselmann  (1976)  shows  that  an  idealized  li¬ 
near  ocean  with  long  response  times  will  redden  the 
spectrum  of  a  white  noise  forcing.  Similarly,  Griffies  and 
Tziperman  (1995)  show  that  high-frequency  stochastic 
atmospheric  forcing  excites  low  frequency  modes  in  the 
ocean.  A  number  of  studies  implicate  intrinsic  nonlinear 
processes,  in  both  the  atmosphere  and  ocean,  as  another 
important  contributor  to  low-frequency  variability  (e.g., 
Holland  1978;  Sheng  and  Hayashi  1990a,b;  Spall  1996; 
Berloff  and  McWilliams  1999;  Dijkstra  and  Molemaker 
1999;  Dewar  2001).  Recently,  Penduff  et  al.  (2011)  dem¬ 
onstrate  that  inter  annual  variance  of  sea  surface  height 
in  ocean  general  circulation  models  driven  by  atmo¬ 
spheric  forcing  fields  that  lack  interannual  variability  is 
comparable  to  the  variance  in  ocean  models  driven  by 
interannually  varying  forcing  fields.  However,  these 
variances  are  comparable  if  and  only  if  the  ocean  models 
are  turbulent  [or,  using  terminology  borrowed  from 
Hecht  and  Hasumi  (2008),  “eddying”].  Penduff  et  al. 
(2011)  and  related  studies  (e.g.,  Taguchi  et  al.  2010) 
therefore  suggest  a  role  for  mesoscale  eddies  in  the 
maintenance  of  oceanic  low-frequency  variability. 

As  pointed  out  by  Ferrari  and  Wunsch  (2010),  the 
theory  of  quasi-two-dimensional  turbulence  is  less  well 
developed  in  the  frequency  domain  than  in  the  wave- 
number  domain.  Numerous  studies  (e.g.,  Fjortoft  1953; 
Batchelor  1953;  Kraichnan  1967,  1971;  Chamey  1971; 
Vallis  2006)  argue  that  nonlinearities  in  two-dimensional 
and  quasigeostrophic  (QG)  turbulence  drive  an  inverse 
cascade  of  energy  toward  larger  length  scales  (lower 
wavenumbers).  Scott  and  Wang  (2005)  document  an 
inverse  cascade  in  the  surface  ocean,  strongly  suggesting 
that  the  ocean  is  in  a  state  of  geostrophic  turbulence. 
Arbic  et  al.  (2012b,  hereinafter  ASFMRS)  demonstrate 
that  nonlinearities  in  an  idealized,  two-layer, /-plane,  QG 
turbulence  model  forced  by  an  imposed,  baroclinically 
unstable,  mean  flow  (e.g.,  Salmon  1978, 1980;  Haidvogel 
and  Held  1980;  Larichev  and  Held  1995)  drive  energy 
toward  longer  time  scales  (lower  frequencies),  alongside 
the  well-known  inverse  cascade  to  lower  wavenumbers. 
ASFMRS  also  show  that  nonlinearities  play  a  larger  role 
than  forcing  and  damping  in  the  maintenance  of  low- 
frequency  variance,  at  least  in  baroclinically  unstable 
/-plane  QG  turbulence.  The  frequency  domain  spectral 
transfers  and  fluxes  used  in  ASFMRS  provide  a  new  tool 


to  measure  the  relative  importance  of  nonlinear  intrinsic 
processes  versus  forcing  in  the  maintenance  of  low- 
frequency  oceanic  variability.  ASFMRS  show  that  non- 
linearities  also  generally  drive  energy  to  longer  time 
scales  in  realistic  eddying  ocean  general  circulation  mo¬ 
dels.  We  note  that  the  term  “forcing”  can  mean  different 
things  in  different  settings.  In  our  idealized  QG  model  the 
forcing  is  an  imposed,  horizontally  homogeneous,  tem¬ 
porally  fixed,  baroclinically  unstable  mean  flow.  In  the 
realistic  models  used  by  Penduff  et  al.  (2011)  and  Taguchi 
et  al.  (2010),  forcing  denotes  the  atmospheric  fields, 
which  vary  in  complex  ways  in  both  space  and  time. 

Here  we  extend  the  frequency  domain  analysis  of 
ASFMRS  to  the  frequency-wavenumber  domain. 
Frequency-wavenumber  analysis  allows  inspection  of 
the  behaviors  at  low  frequencies  versus  low  wave- 
numbers,  a  determination  of  which  wavenumbers  are 
primarily  responsible  for  nonlinearity  at  low  frequen¬ 
cies,  and  a  determination  of  which  frequencies  are  pri¬ 
marily  responsible  for  nonlinearity  at  low  wavenumbers. 
Frequency-wavenumber  spectra  provide  a  test  of  the 
Taylor  (1938)  hypothesis  relating  frequency  spectra  to 
wavenumber  spectra.  Finally,  as  will  be  seen,  the  anal¬ 
ysis  of  the  frequency-wavenumber  spectral  fluxes  pro¬ 
vides  a  constraint  on  the  values  of  the  frequency  domain 
spectral  fluxes  at  the  lowest  frequencies,  which  in  turn 
allows  for  a  closer  comparison  of  the  low-frequency 
fluxes  in  QG  versus  realistic  models  than  was  done 
in  ASFMRS.  Previous  studies  of  oceanic  frequency- 
wavenumber  spectra  include  Wunsch  and  Stammer 
(1995),  Chelton  and  Schlax  (1996),  Farrar  (2008), 
Wunsch  (2009, 2010),  Early  et  al.  (2011),  and  Wortham 
(2013),  among  others.  A  comprehensive  discussion  of 
frequency-wavenumber  spectra  computed  from  models, 
satellite  altimetry  data,  and  in  situ  data  can  be  found 
in  Wortham  (2013).  Hayashi  and  Golder  (1977)  and 
Hayashi  (1982)  discuss  the  “space-time”  (wavenumber- 
frequency)  spectra  of  realistic  atmospheric  models  and 
datasets.  We  are  not  aware  of  any  previous  studies  ex¬ 
amining  frequency-wavenumber  domain  spectral  fluxes 
and  transfers  in  idealized  QG  turbulence  models,  or  in 
an  oceanic  context,  as  we  do  here.  Note  that  our  focus 
here  is  on  isotropic  wavenumber-frequency  spectral 
diagnostics.  In  a  companion  paper,  we  analyze  spectral 
diagnostics  in  the  anisotropic  wavenumber-frequency 
domain  to  better  understand  the  role  of  nonlinearities  in 
westward  eddy  propagation. 

We  examine  output  from  the  same  two-layer  QG 
turbulence  model  utilized  in  ASFMRS.  However,  the 
realistic  eddying  ocean  model  utilized  here  is  the  Hy¬ 
brid  Coordinate  Ocean  Model  (HYCOM;  Chassignet 
et  al.  2007),  in  place  of  the  Naval  Research  Laboratory 
(NRL)  Layered  Ocean  Model  (NLOM;  Hurlburt  and 
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Thompson  1980;  Shriver  et  al.  2007)  used  in  ASFMRS. 
HYCOM  has  more  layers  in  the  vertical  direction  than 
does  NLOM  and  a  number  of  other  more  realistic  features. 

As  in  ASFMRS  and  Arbic  et  al.  (2013),  we  will  also 
compute  spectral  diagnostics  from  a  satellite  altimeter 
product.  The  Archiving,  Validation,  and  Interpretation 
of  Satellite  Oceanographic  (AVISO)  gridded  satellite 
altimeter  product  (Le  Traon  et  al.  1998;  Ducet  et  al.  2000) 
we  use  is  constructed  through  interpolation,  in  space  and 
time,  of  along-track  sea  surface  height  (SSH)  data  from 
two  altimeters  onto  a  Vf  Mercator  grid  at  regular  7-day 
intervals.  Because  interpolation  smooths  the  along-track 
data,  the  spatial  feature  resolution  of  AVISO  maps  is 
coarser  than  the  ~6-km  resolution  along  track  and  is 
closer  to  the  ~300-km  cross-track  spacing  of  the  Ocean 
Topography  Experiment  (TOPEX)/Poseidon  and  Jason 
altimeters  (Chelton  et  al.  2011).  Wortham  (2013)  finds 
that  the  AVISO  product  has  significantly  lower  energy 
than  raw,  along-track  altimetry  data  for  wavelengths 
shorter  than  about  225  km  and  periods  shorter  than  about 
55  days.  We  investigate  the  impact  of  the  smoothing  in¬ 
herent  in  the  construction  of  AVISO  products  through 
comparison  of  spectral  fluxes  computed  from  filtered 
versions  of  the  realistic  model  output  with  fluxes  com¬ 
puted  from  unadulterated  output. 

This  paper  is  organized  as  follows:  In  section  2,  we 
describe  the  models  and  altimeter  product.  In  section  3, 
we  discuss  the  theory  of  frequency-wavenumber  domain 
spectral  fluxes  and  transfers.  Section  4  presents  the  anal¬ 
ysis  technique,  and  section  5  presents  the  results.  A  sum¬ 
mary  and  discussion  is  given  in  section  6. 

2.  Description  of  models  and  gridded  altimeter 
product 

a.  Idealized  quasigeostrophic  turbulence  model 

For  simplicity,  we  utilize  a  highly  idealized  two-layer 
QG  turbulence  model  on  an  /-plane  with  a  flat  bottom, 
forced  by  an  imposed,  horizontally  homogeneous,  ver¬ 
tically  sheared  mean  flow  that  roughly  represents  oceanic 
gyres.  The  model  is  pseudospectral  and  doubly  periodic, 
with  256  X  256  grid  points,  where  each  side  has  di¬ 
mensional  length  Flow  begins  with  a  randomly 

generated  initial  condition,  and  the  model  achieves  sta¬ 
tistical  equilibrium  when  eddy  energy  extraction  from  the 
baroclinically  unstable  mean  flow  is  balanced  by  energy 
dissipation.  Despite  the  numerous  simplifications,  the 
frequency/isotropic  wavenumber  behaviors  in  the  QG 
model  compare  well  in  many  respects  with  the  behaviors 
in  the  realistic  ocean  model  and  in  gridded  satellite  al¬ 
timeter  maps. 

We  follow  the  format  of  Flierl  (1978)  in  our  governing 
equations: 


^  +  +/(i//^,g^)  =  ssd,  and  (1) 

-^+ Fj +7(i/'2,?2)  ^ 

where  ssd  is  small-scale  dissipation.  Here  t  is  time,  q  is 
the  perturbation  potential  vorticity,  and  subscripts  1  and 
2,  respectively,  denote  upper  and  lower  layers.  Pertur¬ 
bation  stre amf unctions  if/i  and  1/^2  are  related  to  q  by 


+ 


i<l'2  ~  '/'l) 

(i+5)l2; 


<72  =  + 


-  l/'2) 
(1+5)L2’ 


(3) 


where  8  =  H1IH2,  the  ratio  of  upper-  to  lower-layer 
depths,  and  the  first  baroclinic  mode  deformation  radius 
Ld  is  defined  by  [g'HiH2lf^{Hi  +  H2)f^,  where  g'  is 
reduced  gravity  and  /o  is  the  Coriolis  parameter.  The 
zonal  (east-west)  and  meridional  (north-south)  spatial 
coordinates  are  x  and  y,  respectively.  The  Jacobian  of 
two  fields  A  and  B  is  given  by  J(A,  B)  =  dA/ 
dx  dB/dy  —  dAldy  dBldx.  The  mean  flow  forcing  terms  are 


and 

dx  dy  dx 

(4) 

dq2  dq^  dil/2 

(5) 

dx  dy  dx  ’ 

where  imposed  time-mean  quantities  are  designated 
with  overbars.  The  zonal  and  meridional  velocities  are 
denoted  by  u  and  u,  respectively.  On  the  /-plane,  the 
direction  of  the  mean  flow  is  arbitrary;  we  take  it  to  be 
zonal.  The  imposed  time-mean  potential  vorticity  gra¬ 
dients  are 

dq^  {u^  —  u^)  dq^  8{u^  —  Iq) 

■^"(1  +  5)L2’ 


The  bottom  Ekman  drag  R2  is  related  to  the  bottom 
boundary  layer  thickness  (^Ekman  =  \/A//o  by 


R 


2 


/p^Ekman 

2^2  ’ 


(7) 


where  A  is  the  vertical  eddy  viscosity.  We  use  an  expo¬ 
nential  cutoff  filter  in  wavenumber  space  (Canute  et  al. 
1988;  LaCasce  1996)  for  ssd: 

filter  =  exp ^ a(^\/ R  -  j  when  + 
filter  =  1.0  when  +  P  ^  , 

(8) 
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where  k  and  I  are  zonal  and  meridional  wavenumbers, 
respectively,  and  is  a  cutoff  wavenumber.  As  in 
LaCasce  (1996),  which  can  be  consulted  for  more  de¬ 
tails,  we  take  a  =  18.4,  M  =  4,  and  ko  =  0.65kN,  where 
kjs/  corresponds  to  the  Nyquist  scale.  At  every  time  step 
the  filter  is  applied  to  qi  and  the  Fourier  transforms 
of  qi  and  ^2- 

Our  QG  equations  are  governed  by  two  nondimensional 
parameters.  We  set  the  stratification  parameter  8  to  0.2, 
which  yields  the  strength  of  the  first  baroclinic  mode 
self-interaction  in  the  two-layer  model  consistent  with 
that  obtained  from  typical  midlatitude  stratification 
profiles  (Flierl  1978;  Fu  and  Flierl  1980).  The  non- 
dimensional  linear  bottom  Ekman  friction  strength, 
Fl  =  R2Ldl(ui  -  ui),  is  set  to  0.4,  which  yields  ampli¬ 
tudes,  horizontal  scales,  and  vertical  structure  of  eddy 
kinetic  energy  (EKE)  that  agree  reasonably  well  with 
midlatitude  observations  (Arbic  and  Flierl  2004;  Arbic 
et  al.  2013).  Figure  1  displays  representative  snapshots  of 
ij/i  and  qi  for  the  QG  model.  The  flow  field  consists  of 
densely  packed  eddies  in  qualitative  similarity  to  a  typi¬ 
cal  snapshot  of  midocean  SSH.  Isolated  coherent  vorti¬ 
ces  (e.g.,  McWilliams  1984)  can  be  seen  in  the  if/i  and 
especially  qi  fields.  The  spectral  diagnostics  are  calcu¬ 
lated  from  5400  snapshots,  each  spaced  one  unit  of 
nondimensional  time  LJijli  —  ui)  apart.  The  correlation 
time  (as  measured  by  the  zero  crossing  in  the  autocor¬ 
relation  of  i/^i)  is  about  16.5  units  of  nondimensional  time, 
suggesting  that  our  sampling  rate  is  adequate. 

b.  Global  eddying  general  circulation  model 

For  details  on  global  HYCOM  the  reader  is  directed 
to  Metzger  et  al.  (2010)  and  references  therein.  A  brief 
summary  is  given  here.  The  HYCOM  simulation  used  in 
this  paper  is  nonassimilative,  has  32  layers  in  the  vertical 
direction,  and  has  a  horizontal  resolution  of  1/12.5°  at 
the  equator.  The  simulation  was  run  interannually  over 
the  period  July  2003-December  2010  using  3-hourly 
Fleet  Numerical  Oceanography  Center  Navy  Opera¬ 
tional  Global  Atmospheric  Prediction  System  (FNMOC 
NOGAPS;  Rosmond  et  al.  2002)  atmospheric  forcing 
with  wind  speeds  scaled  to  be  consistent  with  Quick 
Scatterometer  (QuikSCAT;  e.g.,  Liu  and  Xie  2006)  ob¬ 
servations.  Comparison  to  surface  drifters  and  deep 
current  meters  demonstrates  that  HYCOM  generates 
a  realistic  circulation,  albeit  somewhat  deficient  in  EKE, 
being  only  79%  and  80%  of  the  observations  at  the 
surface  and  abyssal  ocean  (Thoppil  et  al.  2011).  Because 
we  are  limited  by  the  substantial  computational  re¬ 
sources  needed  to  run  global  eddying  HYCOM  and  to 
store  its  output,  we  have  only  5  yr  (2005-09)  of  HYCOM 
SSH  output  for  our  analysis.  Storage  concerns  preclude 
us  from  examining  depth-integrated  energy  budgets  in 


(a) 


Fig.  1.  Representative  snapshots  of  (a)  upper-layer  stream- 
function  ipi  [normalized  by  (wi  -  U2)Ld\  and  (b)  upper-layer  poten¬ 
tial  vorticity  qi  [normalized  by  iui  ~  U2)ILd\  in  idealized  two-layer 
QG  simulation.  The  domain  is  doubly  periodic  and  on  a  side, 

where  is  the  first  baroclinic  mode  Rossby  deformation  radius. 

HY COM,  as  we  do  in  our  QG  analysis.  Instead  we  focus 
on  the  spectra,  spectral  fluxes,  and  spectral  transfers  of 
geostrophic  kinetic  energy  in  HYCOM’s  sea  surface. 
We  choose  to  analyze  a  HYCOM  simulation  forced  by  the 
astronomical  tidal  potential  as  well  as  atmospheric  fields 
(Arbic  et  al.  2010,  2012a;  Shriver  et  al.  2012;  Richman 
et  al.  2012)  so  that  ongoing  frequency-wavenumber 
analyses  of  tidal  and  near-inertial  flows  can  be  easily 
compared  to  the  analyses  of  geostrophic  flows  performed 
here.  A  3-day  low-pass  filter  is  employed  to  remove  high- 
frequency  motions  from  the  hourly  HY  COM  SSH  fields. 
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Fig.  2.  Six  regions  used  to  compute  spectra,  spectral  fluxes,  and  spectral  transfers  of  surface 
ocean  geostrophic  kinetic  energy  from  (a)  1/12.5°  HYCOM  output  and  (b)  V3°  AVISO  gridded 
satellite  altimeter  product.  In  (a)  the  regions  are  highlighted  against  the  0000  UTC  1  January 
2008  snapshot  of  HYCOM  steric  sea  surface  height  (cm).  In  (b)  the  regions  are  highlighted 
against  the  AVISO  20  December  2006  snapshot  of  sea  surface  height  anomaly  (cm). 


SO  that  only  low-frequency  (assumed  geostrophic)  mo¬ 
tions  remain.  We  analyze  low-pass  filtered  results  spaced 
at  1-day  intervals. 

Six  regions  of  interest,  over  which  we  will  compute 
spectral  diagnostics,  are  shown  in  Fig.  2a,  with  a  snapshot  of 
HYCOM  steric  SSH  serving  as  background.  Four  of  the 
regions  experience  intense  boundary  current  and  mesoscale 
eddy  activity — the  Kuroshio  in  the  western  North  Pacific, 
the  Gulf  Stream  in  the  western  North  Atlantic,  the  Agulhas 
off  the  southern  tip  of  Africa,  and  the  Malvinas  in  the 
western  South  Atlantic.  Two  additional,  relatively  quiescent 
regions  are  in  the  southeastern  Pacific — a  midlatitude  re¬ 
gion  and  a  high-latitude  region  that  touches  upon  the 
Antarctic  Circumpolar  Current  in  its  southernmost  por¬ 
tions.  Because  of  the  removal  of  HYCOM  land  points  as¬ 
sociated  with  Bermuda,  our  Gulf  Stream  region  is  slightly 
narrower  in  east-west  extent  in  HYCOM  than  in  AVISO. 

c.  Gridded  satellite  altimeter  product 

We  use  840  weekly  snapshots  of  the  two-satellite 
AVISO  V3°  Mercator  grid  reference  product  (Le  Traon 


et  al.  1998;  Ducet  et  al.  2000),  beginning  14  October  1992 
and  ending  12  November  2008.  Figure  2b  displays  the 
six  regions  over  which  we  compute  AVISO  spectral 
diagnostics,  shown  against  an  example  AVISO  SSH 
anomaly  snapshot.  The  regions  correspond  very  closely 
to  those  utilized  in  the  HYCOM  analysis,  with  the  greatest 
deviation  being  in  the  Gulf  Stream  as  discussed  above. 

d.  Normalization  of  wavenumbers  and  frequencies 

We  end  this  section  with  a  brief  discussion  of  the 
X  axes  used  in  spectral  plots  throughout  the  paper.  In 
the  QG  plots,  wavenumbers  are  nondimensionalized 
by  1/Lj  and  frequencies  are  nondimensionalized  by 
f//Lj,  where  1/  =  -  W2  is  the  imposed  mean  shear.  To 

connect  the  HYCOM  and  AVISO  results  with  the  QG 
results,  we  will  display,  on  the  HYCOM  and  AVISO 
plots,  estimated  oceanic  values  of  the  l/L^  wave- 
number  and  UILd  frequency.  For  each  of  the  six  regions 
shown  in  Fig.  2,  we  choose  the  value  at  the  Chelton  et  al. 
(1998)  atlas  grid  point  lying  closest  to  the  center  of  the 
region.  The  value  of  U  is  computed  as  the  area-weighted 


August  2014 


ARBIC  ET  AL. 


2055 


Table  1.  Values  of  L^,  t/,  and  L^IU  (see  text)  for  the  six  regions 
shown  in  Fig.  2, 


Region 

Ld  (km) 

U  (cm  s  ^) 

LdIU 

(days) 

Midlatitude  southeast  Pacific 

36 

2.4 

17 

High-latitude  southeast 

Pacific 

21 

4.6 

5.2 

Agulhas 

30 

13 

2.5 

Malvinas 

25 

9.4 

3.0 

Gulf  Stream 

34 

9.6 

4.1 

Kuroshio 

32 

8.8 

4.2 

root-mean-square  of  the  time-mean  surface  velocity  over 
the  region: 


(9) 


where  dA  represents  an  element  of  area,  and  Wmean  and 
i^mean  are  respectively  the  zonal  and  meridional  surface 
velocities  computed  from  the  mean  dynamic  topogra¬ 
phies  of  Niiler  et  al.  (2003)  and  Maximenko  et  al.  (2009). 
We  assume  that  the  surface  velocity  is  a  reasonable 
proxy  for  the  shear,  that  is,  that  velocities  in  the  sub¬ 
surface  ocean  are  substantially  weaker  than  the  surface 
flow.  The  values  of  L^,  U,  and  LJU  for  the  six  regions 
are  given  in  Table  1. 

3.  Theory  of  frequency-wavenumber  spectral 
fluxes  and  transfers 

The  theory  of  spectral  transfers  T{K,  w)  and  spectral 
fluxes  Il{K,  (o)  in  the  two-layer  QG  model,  where 
k=VWTp  is  isotropic  wavenumber  and  co  denotes 
frequency,  is  summarized  here.  We  let  A  denote  the 
discrete  Fourier  transform  of  an  arbitrary  held  A.  For 
example, 

'Pi{x,y,  0  =  S  II  Mk,  I,  .  (10) 

k  I  (O 

Thus,  the  frequency-wavenumber  spectrum  of  upper- 
layer  kinetic  energy  density  |Vr/^ip/2  is 

EKEj (k,  I,  (o)  =  ^K^^{k,  I,  (t))i^(L  l,(o),  (11) 

where  the  superscript  *  represents  a  complex  conjugate. 

As  discussed  in  many  sources  (e.g.,  Salmon  1978, 1980; 
Hua  and  Haidvogel  1986;  Larichev  and  Held  1995; 
LaCasce  1996),  spectral  transfers  identify  energy  sources 
and  sinks  in  Fourier  space,  which  could  mean  wavenumber, 
frequency,  or  frequency-wavenumber  space  here.  Positive 


(negative)  transfer  values  indicate  an  energy  source 
(sink).  To  develop  a  depth-averaged  spectral  transfer 
budget,  we  multiply  the  Fourier  transform  of  the  upper- 
layer  equation  [Eq.  (1)]  by  — (/c,  /,  (o)l{l  +  5),  mul¬ 
tiply  the  Fourier  transform  of  the  lower-layer  equation 
[Eq.  (2)]  by  /,  (o)l{l  +  5),  and  add  the  results, 

yielding 

^KE,2(^’  ^APe(^’ 

^forcmg(L/,^)+r«etion(L/,«)  =  0.  (12) 


Here  the  spectral  transfers  of  upper-  and  lower-layer 
kinetic  energy  are 


^KE  1  ^ 


1  +  8 


/,  (o)J(iljp  V^i/^^)(A:,  /,  (o) 


and 

(13) 


1+5 


/,  co)J(ip2,  /,  (o) 


(14) 


respectively,  where  Re  denotes  the  real  part  of  a  complex 
number.^  The  spectral  transfers  of  available  potential 
energy  (APE),  mean  flow  forcing,  and  bottom  Ekman 
friction  are 


^ape(^’  (o)  =  Re 


{i+sfq 

I,  <w) 


(15) 


+rcing(^’^>")  =  Re 


1+5 
1 

1  +  5 


^i{k,  /,  (o)T^{k,  /,  (o) 


,  and 


(16) 


fiction  k  «) , 


(17) 


^  For  simplicity,  because  stratification  is  not  horizontally  uniform 
in  HYCOM  and  AVISO  output,  the  factor  5/(1  +  8)  is  omitted 
from  the  HYCOM  and  AVISO  spectral  kinetic  energy  flux  and 
transfer  computations  in  this  paper. 
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respectively.  The  transfer  terms  associated  with  ssd,  not 
listed  in  Eq.  (12),  are  inferred  as  a  small  residual.  Note 
that  because  Eq.  (12)  is  a  depth-averaged  budget  there 
are  no  interlayer  transfers. 

Spectral  fluxes  are  defined  as  integrals  of  the  spectral 
transfers,  namely, 

nterm  =  f  f  f  Lerm  ")  ’ 

(18) 


where  the  subscript  “term”  denotes  KE,  1;  KE,  2;  APE; 
forcing;  or  friction  as  appropriate.^  Integration  over  all 
frequencies  co  yields  spectral  fluxes  in  isotropic  wave- 
number  K  space 


nte™(^) 


^term(^’  ’ 


(19) 


while  integration  over  all  wavenumbers  yields  spectral 
fluxes  in  frequency  space 


nterm(fi) 


^term 

j  J  Jw>n 


4.  Analysis  technique 

Spectral  analysis  of  fields  that  are  not  periodic  re¬ 
quires  care  (Priestley  1981).  None  of  the  three  output 
fields  analyzed  here — from  the  idealized  QG  model, 
realistic  model,  or  altimeter  product — are  periodic  in 
time,  nor  are  fields  excised  from  the  realistic  model  or 
altimeter  product  periodic  in  space.  For  the  AVISO  and 
3-day  low-pass  filtered  HYCOM  output,  we  first  convert 
the  SSH  fields  ry  to  surface  geostrophic  streamfunctions, 
denoted  by  ij/i  in  analogy  to  the  QG  model,  via  ij/i  =  grj/f, 
where  g  =  9.8  m  s~^  is  the  gravitational  acceleration,  and 
/  is  the  Coriolis  parameter.  We  remove  the  temporal 
trend  and  mean  of  if/i  at  each  spatial  grid  point.  Fol¬ 
lowing  this,  spatial  trends  and  means  are  removed  from 
every  snapshot  of  ij/i  with  a  two-dimensional  least 
squares  fit.  To  enforce  periodicity  in  space,  each  if/i 
snapshot  is  multiplied  by  a  two-dimensional  Tukey 
window  having  a  ratio  of  taper  to  constant  sections  equal 
to  0.2  in  each  horizontal  direction.  To  enforce  period¬ 
icity  in  time,  the  ij/i  time  series  at  each  spatial  grid  point 


is  multiplied  by  a  Tukey  window,  also  having  a  value  of 
0.2  for  the  ratio  of  taper  to  constant  sections.^  Because 
the  QG  model  is  doubly  periodic,  and  because  if/i  has 
a  zero  time  mean  at  each  grid  point,  we  do  not  employ 
any  spatial  detrending  and  mean  removal,  spatial  win¬ 
dowing,  or  temporal  detrending  and  mean  removal  in 
our  QG  analysis.  However,  the  QG  analysis  does  em¬ 
ploy  a  temporal  Tukey  window  to  enforce  periodicity  in 
time.  Terms  involving  spatial  derivatives  are  computed 
using  forward  and  inverse  discrete  Fourier  transforms  in 
wavenumber  space.  For  example,  is  the  inverse 
discrete  Fourier  transform  of  +  For 

simplicity,  in  the  HYCOM  and  AVISO  analyses  we 
neglect  the  small  convergence  of  meridians  within  our 
regions.  We  use  the  native  HYCOM  and  AVISO  grids, 
rather  than  interpolating  to  regularly  spaced  x-y  grids, 
neglecting  the  order  10%  variations  in  latitudinal  spac¬ 
ing  on  Mercator  grids. 

Multiplication  by  a  Tukey  window,  after  the  removal 
of  a  trend  and  mean,  leaves  behind  fields  having  nonzero 
means.  Therefore,  some  nonzero  spectral  energy  is 
found  in  the  zero  frequency  (wavenumber).  As  is  com¬ 
monly  done,  we  will  not  include  the  zero  frequency  and 
wavenumber  results  in  our  plots  of  spectra.  However,  as 
will  be  seen  shortly,  the  zero-frequency  results  are  cru¬ 
cial  for  spectral  fluxes  and  transfers.  To  display  the  zero- 
frequency  (wavenumber)  flux  and  transfer  results  on  the 
semilog  plots  we  favor  throughout  this  paper,  we  move 
the  zero-frequency  (wavenumber)  results  to  one-half  of 
the  lowest  resolved  nonzero  frequency  (wavenumber). 

As  described  below,  consideration  of  a  joint  frequency- 
wavenumber  spectral  flux  provides  insights  into  the  low- 
frequency  behavior  of  nKE,i(<^)*  For  geostrophic  flows 
with  periodic  boundary  conditions,  nKE,i(^)  is  zero  at 
the  lowest  wavenumber.  Because  HkeaC^^  is  a  dou¬ 
ble  integral  over  wavenumber  and  frequency  of 
Fke,i(^5  because  the  order  of  integration  should 

not  matter,  nKE,i(^)  niust  be  zero  at  the  lowest  fre¬ 
quency.  For  all  of  the  nKE,i(^5  curves  in  this  paper, 
the  absolute  value  at  the  lowest  frequency-lowest 
wavenumber  pair  is  2%  or  less,  often  much  less,  of  the 
largest  absolute  value  occurring  throughout  (K,  co) 
space.  Because  nKE,i(^)  rnust  be  zero  at  the  lowest 
frequency,  the  record  duration  impacts  estimated  low- 
frequency  nKE,i(<^)  values.  In  Fig.  3a,  we  display  the 


^In  retrospect,  we  realize  that  the  factor  of  6/(1  +  8)  was 
sometimes  inadvertently  omitted  from  ASFMRS  plots  in  which  the 
QG  nKE,i(^)  curve  was  shown  as  a  stand-alone  result — that  is,  not 
as  part  of  a  depth-averaged  budget.  Because  this  omission  was  not 
made  in  any  of  the  energy  budget  plots,  it  does  not  affect  any  of  the 
conclusions  of  ASFMRS. 


^Note  that  most  of  the  analysis  in  ASFMRS  and  Arbic  et  al. 
(2013)  utilized  nine  overlapping  Hanning  windows  in  space  (and 
time,  in  the  case  of  ASFMRS),  instead  of  the  Tukey  windows  used 
here.  As  discussed  in  Arbic  et  al.  (2013)  overlapping  Hanning 
windows  and  Tukey  windows  yield  quantitatively  similar  flux  es¬ 
timates  in  high-resolution  models  and  qualitatively  similar  AVISO 
fluxes. 
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(a)  Impact  of  simulation  duration  on  QG  ^(co)  values 


_ 
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10-^  10-^  10-^  10°  10^ 
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Fig.  3.  (a)  Spectral  kinetic  energy  fluxes  nominal  QG 

simulation  (Long;  black)  and  shorter  subset  of  the  QG  simulation 
(Short;  red).  In  the  axis  labels  of  (a)  and  of  other  figures  made  from 
the  QG  simulation,  “U”  denotes  the  imposed  mean  shear  ui-U2. 
(b),(c)  Spectral  kinetic  energy  fluxes  and  transfers  nKE,i(^) 
Lke,i(^)^  respectively,  in  the  Agulhas  region  of  HYCOM  for  three 
analysis  procedures  given  in  legend  (see  text). 


nKE,i(^)  curves  from  our  nominal  (5400  snapshot)  QG 
simulation  and  from  a  600-snapshot  subset  of  this 
simulation  (having  the  same  temporal  sampling  in¬ 
terval).  Because  nKE,i(<^)  goes  to  near-zero  values 
at  the  low-frequency  ends  of  both  curves,  and  because 
the  low-frequency  end  points  differ  between  the  two 
curves,  the  nKE,i(^)  curves  separate  slightly  at  low 
frequencies. 

The  importance  of  temporal  detrending  is  illustrated 
in  Figs.  3b  and  3c.  The  HY COM  Agulhas  region  nKE,i(^) 
curve  resulting  from  our  preferred  analysis  procedure 
(Detrend) — in  which  a  temporal  mean  and  linear  trend 
is  removed  from  each  spatial  grid  point — is  contrasted 
with  results  from  two  other  procedures  in  Fig.  3b.  The 
High-pass  curve  uses  a  high-pass  to  remove  frequencies 
lower  than  twice  the  lowest  resolved  frequency.  In  the 
Neither  curve,  neither  a  high-pass  or  detrend/mean  re¬ 
moval  is  employed  on  the  time  series  at  each  model  grid 
point.  While  the  Neither  curve  displays  some  positive 
nKE,i(<^)  values  at  low  frequencies,  the  High-pass  and 
Detrend  curves  display  negative  values  (implying  a  cas¬ 
cade  to  lower  frequencies)  throughout  the  entire  plotted 


frequency  range.  The  nonlinear  spectral  kinetic  energy 
fluxes  nKE,i(<^)  ^^c  strongly  affected  at  low  frequencies 
if  temporal  means  and  trends  are  not  removed,  consis¬ 
tent  with  the  well-known  effects  on  the  spectra  them¬ 
selves  (e.g.,  Priestley  1981).  Because  detrending  is 
standard  practice  in  time  series  analysis,  we  believe  that 
the  Detrend  and  High-pass  curves  in  Fig.  3b  are  more 
reliable  than  the  Neither  curve.  Because  our  intuition 
was  developed  in  the  QG  case,  for  which  these  consid¬ 
erations  are  less  important  as  described  earlier,  we  did  not 
recognize  the  need  for  temporal  detrending  in  our 
ASFMRS  analyses,  nor  did  we  assign  zero-frequency  flux 
and  transfer  values  to  displayable  nonzero  frequencies. 
As  a  result  the  realistic  model  nKE,i(<^)  curves  in 
ASFMRS  look  like  the  Neither  curve  in  Fig.  3b,  but 
without  the  lowest  frequency  result — in  other  words,  the 
ASFMRS  realistic  model  nKE,i(<^)  results  show  artificial 
positive  values  at  low  frequencies,  without  a  drop  to  near¬ 
zero  values  at  the  left-hand  end  point.  The  impact  of  these 
important  signal  processing  concerns  on  the  HYCOM 
Agulhas  rKE,i(<^)  curves  is  shown  in  Fig.  3c.  At  the  lowest 
frequencies,  the  Neither  curve  has  an  unphysically  sharp 
jump  due  to  substantial  energy  in  the  (nominally)  zero 
frequency.  In  the  remainder  of  this  paper,  we  employ  the 
Detrend  procedure. 


5.  Results 

a.  Spectra 

We  now  discuss  frequency-wavenumber  spectra 
EKEi  of  ocean  surface  kinetic  energy  [Eq.  (11)].  The 
ratio  of  HYCOM  to  AVISO  total  kinetic  energy  (in¬ 
tegrated  over  all  frequencies  and  wavenumbers)  is  of 
order  one — 0.85  in  the  Kuroshio  region  and  between  1.1 
and  1.3  in  the  other  five  regions.  Frequency-wavenumber 
spectra  for  the  QG  model  and  for  HYCOM  and  AVISO 
results  in  the  Agulhas  region,  which  is  reasonably  rep¬ 
resentative,  are  shown  in  Fig.  4.  The  spectra  display 
several  visual  similarities,  including  a  tendency  for  high- 
wavenumber  energy  to  be  nearly  independent  of  fre¬ 
quency,  a  tendency  for  a  flat  frequency  spectrum  at  low 
frequencies,  and  a  tendency  for  lower-wavenumber  en¬ 
ergy  to  be  associated  with  lower  frequencies.  This  latter 
feature  is  seen  as  a  diagonal  “cut”  in  the  right-hand  side 
of  each  subplot  (implying  small  low-wavenumber  energy 
at  high  frequencies).  The  spectra  are  visually  “wider”  in 
frequency  space  than  in  wavenumber  space.  Substantial 
values  of  EKEi  are  seen  over  a  decade  or  less  in  wave- 
number,  but  over  two  or  more  decades  in  frequency. 
Because  EKEi  rolls  off  at  the  lowest  wavenumbers, 
analyses  in  larger  domains  would  not  see  a  continuation 
of  substantial  EKEi  values  out  to  yet  lower  wavenumbers. 
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(a)  QG 


0 


CO  (rai/day) 
(c)  AVISO  Agulhas 


Fig.  4.  Wavenumber-frequency  {K-oj)  spectrum  EKEi(/^,  oj)  of 
|Vj//ip/2  (kinetic  energy)  from  the  (a)  upper  layer  of  the  idealized 
two-layer  QG  simulation  and  from  Agulhas  region  of  ocean  surface 
in  (b)  HYCOM  and  (c)  AVISO.  Units  of  spectra  are  logio  of 
0.5(wi  -  W2)^/([rad/L^i]{rad/[L£//(wi  -  W2)]})  in  (a)  and  logio  of 
(cms“^)^/[(radday“^)(radkm“^)]  in  (b)  and  (c).  Dashed  horizon¬ 
tal  and  vertical  lines  respectively  represent  the  deformation 
wavenumber  1  rad/L^  and  1  rad/(L^/t/),  where  U  is  the  mean  shear 
(see  text).  The  x  and  y  axes,  and  the  colorbar,  in  the  HYCOM 
subplot  (b)  are  identical  to  those  in  the  AVISO  subplot  (c).  The  x 
and  y  axes,  and  the  colorbar,  in  the  QG  subplot  (a)  cover  the 
same  number  of  decades  as  do  those  in  the  HYCOM  and  AVISO 
subplots. 


In  contrast,  because  of  the  flatness  of  EKEi  at  low  fre¬ 
quencies,  the  record  duration  determines  the  range  of 
frequencies  over  which  substantial  EKEi  values  are 
seen.  The  QG  record  is  the  longest  of  the  three  shown  in 
Fig.  4 — 86  yr  if  one  takes  typical  values  of  L ^  =  50  km 
and  U  =  lOcms”^ — hence,  it  has  the  largest  width  in 
frequency.  The  asymmetry  between  spectral  widths  in 
wavenumber  versus  frequency  will  be  seen  throughout 
this  paper. 

To  further  investigate  the  association  of  low-frequency 
and  low-wavenumber  energy,  we  compute  the  energy- 
containing  wavenumbers  Ke{(o),  and  energy-containing 
frequencies  (Oe{K),  defined  as  first  moments  (centroids)  of 
the  spectrum,  namely. 


KJo,)  = 


JJ  \/W^'EKEyk,l,m)dkdl 

'^^BKEyk,l,io)dkdl 


(21) 


(for  a  fixed  m)  where  EKEi(/c,  /,  m)  is  the  kinetic  energy 
spectrum  [Eq.  (11)],  and 

mEKE^(i^,  (o)  dco 

coe{K)  =  ^ -  (22) 

EKE^(i^,m)  dco 

for  a  fixed  isotropic  wavenumber  K.  A  general  associa¬ 
tion  between  low-wavenumber  and  low-frequency  en¬ 
ergy  is  found  in  plots  of  (Oe(K)  and  Ke((o)  for  the  models 
and  altimeter  maps  (Fig.  5).  The  correspondence  is 
not  simple,  and  the  plots  are  noisy,  especially  near 
the  Nyquist  wavenumbers  and  frequencies.  Note  that 
a  general  association  of  low-frequency  energy  with  low 
wavenumbers  does  not  preclude  some  association  of 
low-frequency  energy  with  high  wavenumbers.  As  an 
example,  Wortham  (2013)  points  out  that  the  18-yr 
mean  of  AVISO  SSH  contains  considerable  small- 
scale  (high  wavenumber)  structure. 

The  frequency-wavenumber  spectra  shed  light  on  the 
limitations  of  the  Taylor  (1938)  hypothesis  with  respect 
to  geostrophic  turbulence.  The  Taylor  hypothesis  states 
that  when  turbulent  velocities  are  small  compared  to 
a  mean  flow  V,  the  wavenumber  and  frequency  domains 
are  related  simply,  via  co  =  Vk.  Variants  of  Taylor’s 
hypothesis  have  taken  U  to  be  a  root-mean-square 
eddy  speed  (Tennekes  1975)  or  an  eddy  propagation 
speed  (Ferrari  and  Wunsch  2010).  If  the  Taylor  hy¬ 
pothesis  held  exactly,  all  of  the  energy  in  the  frequency- 
wavenumber  spectra  would  be  along  a  line  with  slope 
determined  by  V.  In  this  ideal  case,  integrating  in  either 
wavenumber  or  frequency  would  yield  identical  results, 
because  the  same  point  on  the  line  would  be  picked  up  in 
either  case.  In  fact  substantial  energy  does  lie  along 
nearly  straight  lines  in  frequency-wavenumber  spectra 
(e.g.,  Wunsch  2009,  2010;  Early  et  al.  2011;  Wortham 
2013).  However,  for  every  wavenumber  there  is  a  spec¬ 
trum  of  energy  spread  among  all  frequencies,  and  for 
every  frequency  there  is  a  spectrum  of  energy  spread 
among  all  wavenumbers  (Fig.  4).  The  asymmetries  be¬ 
tween  spectra  in  the  two  axis  directions  imply  that  the 
frequency  and  wavenumber  spectra,  as  well  as  their 
slopes,  are  not  easily  relatable.  In  ASFMRS  we  noted 
further  that  the  difference  in  the  slopes  of  stream- 
function  variance  (not  shown  in  the  present  paper)  and 
geostrophic  velocity  variance  (kinetic  energy)  in  the 
frequency  domain  is  —0.8,  not  —2  as  it  is  in  the  wave- 
number  domain  and  as  it  therefore  would  be  in 
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Fig.  5.  (left)  Centroids  (jOe(K)  of  kinetic  energy  spectrum  vs  wavenumber  K.  (right)  Centroids  Ke{co)  of  kinetic 
energy  spectrum  vs  frequency  oj.  Results  shown  for  (a),(b)  upper  layer  of  idealized  two-layer  QG  model,  (c),(d) 
ocean  surface  in  HYCOM,  and  (e),(f)  ocean  surface  in  AVISO.  In  (c),(d),(e),  and  (f),  colors  represent  regions  shown 
in  Fig.  2.  Cyan  is  midlatitude  southeast  Pacific,  green  is  high-latitude  southeast  Pacific,  red  is  Agulhas,  blue  is 
Malvinas,  magenta  is  Gulf  Stream,  and  black  is  Kuroshio.  Dashed  vertical  lines  in  (a),(c),  and  (e)  represent  the 
deformation  wavenumber  1  rad/Lj  (see  text).  Dashed  lines  in  (b),(d),  and  (f)  represent  1  rad/(L^/t/),  where  U  is  the 
mean  shear  (see  text). 
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the  frequency  domain  if  the  Taylor  hypothesis  held 
precisely. 

Wavenumber  spectra,  integrated  over  all  frequencies, 
are  shown  in  Fig.  6  for  the  QG  simulation  and  for 
Agulhas  HYCOM  and  AVISO  output.  Slopes  of  -Vs 
and  —3,  the  predicted  values  for  energy  and  enstrophy 
cascades,  respectively  (e.g.,  Vallis  2006),  are  indicated  by 
dashed  slanted  green  lines.  The  QG  spectrum  (Fig.  6a) 
drops  off  steeply  at  the  highest  wavenumbers  due  to  the 
action  of  the  wavenumber  filter.  A  long  band,  covering 
wavenumbers  slightly  lower  than  the  1/L^  wavenumber 
up  to  higher  wavenumbers  just  before  the  action  of  the 
wavenumber  filter,  sees  a  slope  somewhat  steeper  than 
-3.  A  narrow  band  of  wavenumbers,  lying  at  scales 
somewhat  larger  than  L^,  exhibits  a  shallower  slope  and 
contains  the  largest  EKEi(i^)  values.  At  the  lowest 
wavenumbers  the  spectrum  rolls  off  quickly  to  relatively 
low  values.  The  HYCOM  and  AVISO  wavenumber 
spectra  (Fig.  6b)  resemble  the  QG  spectrum  in  many 
respects.  For  instance,  the  scale  again  lies  near  the 
low-wavenumber  end  of  the  steep-slope  band,  and 
the  spectra  roll  off  rapidly  at  the  lowest  wavenumbers. 
The  HYCOM  and  AVISO  spectra  lie  close  together 
over  much  of  the  range  displayed,  but  at  higher  wave- 
numbers  the  AVISO  spectrum  rolls  off  rapidly  due  to 
the  filtering  inherent  in  the  construction  of  AVISO.  At 
high  wavenumbers,  the  HYCOM  spectrum  is  affected 
by  scale-selective  dissipation  acting  near  the  model  grid 
spacing,  as  well  as  by  the  3-day  low-pass  filter;  in  general, 
high  frequencies  are  associated  with  high  wavenumbers. 
In  this  location  the  HYCOM  wavenumber  spectral 
slope  lies  reasonably  close  to  the  predicted  -3  value  for 
QG  “interior  mode”  turbulence.  Richman  et  al.  (2012), 
Sasaki  and  Klein  (2012),  and  references  therein  pro¬ 
vide  detailed  discussions  of  wavenumber  spectra  in 
realistic  models. 

Frequency  spectra,  integrated  over  all  wavenumbers, 
are  displayed  in  Fig.  7a  for  the  QG  simulation  and  in 
Fig.  7b  for  Agulhas  HYCOM  and  AVISO  output.  A 
slope  of  -2,  the  value  predicted  by  energy  cascade  ar¬ 
guments  applied  to  frequency  space  (Tennekes  and 
Lumley  1972),  is  indicated  by  dashed  slanted  green  lines. 
All  three  spectra  are  flat  at  lower  frequencies  and 
steeply  sloped  at  higher  frequencies,  with  the  LJU  time 
scale  lying  in  the  middle  of  the  steep  slope  band.  This 
basic  shape  is  found  in  frequency  spectra  computed  from 
current  meter  records  (e.g.,  Richman  et  al.  1977;  Ferrari 
and  Wunsch  2010;  Wortham  2013)  and  in  previous 
studies  using  altimeter  products  (e.g.,  Wunsch  2009, 
2010;  Wortham  2013).  The  HYCOM  and  AVISO 
spectra  are  close  together  over  much  of  the  range  shown, 
with  the  largest  deviations  seen  at  high  frequency  where 
the  filtering  inherent  in  the  construction  of  AVISO 


Fig.  6.  Wavenumber  domain  kinetic  energy  spectra  EKEi(K), 
integrated  over  all  cj.  (a)  Upper-layer  kinetic  energy  spectra  in 
idealized  two-layer  QG  simulation.  Dashed  vertical  line  denotes 
the  deformation  wavenumber  1  rad/L^.  (b)  Ocean  surface  kinetic 
energy  spectra  in  HYCOM  and  AVISO  Agulhas  region  output. 
Dashed  vertical  line  denotes  the  deformation  wavenumber  1  rad/L^ 
for  the  Agulhas  region  (see  text).  Dashed  slanted  green  lines 
indicate  slopes  of  -  Vs  and  -3. 

yields  a  steep  rolloff.  HY COM  frequency  spectra  slopes 
over  the  14-68-day  band  (indicated  by  red  vertical  lines 
in  Fig.  7b),  which  lies  within  the  steep  slope  part  of  the 
spectrum,  are  between  —1.6  and  —2.1  over  the  six  re¬ 
gions.  This  range  is  close  to  the  -1.6  to  -2.2  range  seen 
over  the  same  six  regions  in  the  NLOM  realistic  simu¬ 
lations  analyzed  by  ASFMRS  and  to  the  -1.8  slope  in 
the  steeper  part  of  the  QG  spectrum  (ASFMRS),  as  well 
as  to  the  -2  slopes  seen  in  Wortham  (2013) ’s  analyses  of 
models,  along-track  satellite  altimeter  data,  and  in  situ 
observations. 

b.  Spectral  fluxes 

We  discuss  spectral  fluxes  before  spectral  transfers 
because  fluxes  are  smoother  and  hence  provide  an  easier 
platform  for  comparison  of  the  QG,  HYCOM,  and 
AVISO  results.  The  inverse  kinetic  energy  cascade  to¬ 
ward  lower  wavenumbers,  indicated  by  negative  values 
along  the  middle  of  the  wavenumber  axes,  is  readily 
apparent  in  the  QG,  Agulhas  HYCOM,  and  Agulhas 
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(a)  QG  kinetic  energy  spectrum  vs  co,  integrated  over  all  K 


Fig.  7.  Frequency  domain  kinetic  energy  spectra  EKEi(ai)  in¬ 
tegrated  over  all  K.  (a)  Upper-layer  kinetic  energy  spectra  in  ide¬ 
alized  two-layer  QG  simulation.  Dashed  vertical  line  denotes  the 
frequency  lYSLdl(LJU),  where  U  is  the  imposed  mean  shear 
ui  —  U2.  (b)  Ocean  surface  kinetic  energy  spectra  in  HYCOM  and 
AVISO  Agulhas  region  output.  Dashed  vertical  black  line  denotes 
the  frequency  1  rad/(L^/t/),  where  is  the  deformation  radius, 
and  U  is  the  roughly  estimated  mean  shear  for  the  Agulhas  region 
(see  text).  Red  vertical  lines  in  (b)  indicate  periods  of  14  and  68 
days.  Slopes  of  the  HYCOM  frequency  spectra  computed  over  this 
band  are  discussed  in  the  text.  Dashed  slanted  green  lines  indicate 
slope  of  -2. 


AVISO  nKE,i(^’  niaps  shown  in  Fig.  8.  The  positive 
AVISO  flux  values  at  higher  wavenumbers  are  discussed 
extensively  in  Arbic  et  al.  (2013).  The  negative  IIke,! 
values  along  the  frequency  axes  of  Fig.  8  imply  a  non¬ 
linear  kinetic  energy  cascade  toward  lower  frequencies 
alongside  the  inverse  cascade  to  lower  wavenumbers.  As 
with  the  spectra,  the  spectral  fluxes  are  wider  along 
frequency  axes  than  along  wavenumber  axes.  As  will  be 
seen  shortly,  with  the  exception  of  the  midlatitude 
southeast  Pacific,  HYCOM  IIke,!  results  in  other  re¬ 
gions  are  similar  to  the  HYCOM  Agulhas  IIke,!  results, 
while  AVISO  results  are  less  consistent  between  regions 
than  the  HYCOM  results  are. 

The  depth-averaged  budget  of  QG  spectral  fluxes  is 
shown  versus  wavenumber  in  Fig.  9a  and  versus  fre¬ 
quency  in  Fig.  9b.  In  both  subplots  the  lower-layer 
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Fig.  8.  Nonlinear  spectral  kinetic  energy  flux  nKE,i(^^  i^ 
(a)  upper  layer  of  idealized  two-layer  QG  simulation  and  from 
Agulhas  region  of  (b)  HYCOM  and  (c)  AVISO  ocean  surface 
output.  The  units  are  [(wi  -  U2flLd\  in  (a)  and  (nW/kg)  in  (b)  and 
(c).  Dashed  horizontal  and  vertical  lines  are  as  in  Fig.  4.  The  x  and  y 
axes  are  also  as  in  Fig.  4. 


kinetic  energy  and  residual  terms  are  much  smaller  than 
the  other  terms.  Figure  9b  indicates  that  the  kinetic 
energy  cascade  toward  lower  frequencies  is  approxi¬ 
mately  balanced  by  an  APE  cascade  toward  higher 
frequencies,  similar  to  the  wavenumber-space  balance 
between  inverse  kinetic  energy  cascades  and  forward 
APE  cascades  (Salmon  1980)  seen  in  Fig.  9a.  The  posi¬ 
tive  (negative)  value  of  Ftforcmg  (nMction)  at  the  lowest 
wavenumber  and  lowest  frequency  implies  that  the 
forcing  (friction)  adds  energy  to  (removes  energy  from) 
the  system.  Significant  IIke,!  values  cover  about  one 
decade  in  wavenumber  and  two  decades  in  frequency. 

The  nonlinear  surface  ocean  spectral  kinetic  energy 
fluxes  nKE,i(^)  and  nKE,i(<^)  are  shown  across  all  six 
regions  in  Figs.  10  and  11,  respectively,  for  un¬ 
adulterated  HYCOM  output,  HYCOM  output  that  has 
undergone  further  filtering,  and  AVISO  output.  Spatial 
filtering,  temporal  filtering,  and  combined  spatial  and 
temporal  filtering  are  performed  on  HYCOM  output  to 
ascertain  the  effects  of  the  filtering  inherent  in  the 
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(a)  QG  n(K)  budget 


(b)  QG  n(co)  budget 


Fig.  9.  (a)  Depth-averaged  budget  of  wavenumber  domain 
spectral  fluxes  n(A^)  in  QG  model.  Dashed  vertical  line  represents 
the  deformation  wavenumber  1  rad/L  j.  (b)  Depth-averaged  budget 
of  frequency  domain  spectral  fluxes  n(a;)  in  QG  model.  Dashed 
vertical  line  represents  the  frequency  1  rad/(Lj/t/),  where  U  is  the 
imposed  mean  shear  ui  -U2.  Budget  terms  identified  by  legend 
shown  in  (a). 

creation  of  AVISO  products  on  spectral  flux  estimates. 
The  filtering  is  applied  to  HYCOM  output  after  the 
model  has  been  run.  We  utilize  Blackman  filters,  with 
filter  scales  inspired  by  the  spatial  and  temporal  scales  in 
the  AVISO  product.  Because  the  AVISO  product  is  put 
out  every  7  days,  on  a  V3°  Mercator  grid,  we  choose 
a  temporal  Blackman  filter  that  goes  to  zero  7  days  from 
the  central  time  and  a  spatial  Blackman  filter  that  goes 
to  zero  37  km  (V3°  latitude)  from  the  central  grid  point. 

As  in  previous  studies  (Scott  and  Wang  2005;  Scott 
and  Arbic  2007;  Schlosser  and  Eden  2007;  Tulloch  et  al. 
2011;  ASFMRS;  Arbic  et  al.  2013),  negative  nKE,i(^) 
values,  representing  inverse  kinetic  energy  cascades  to 
lower  wavenumbers,  are  ubiquitous  (seen  in  every  sub¬ 
plot  of  Fig.  10).  As  noted  in  Scott  and  Wang  (2005),  the 
deepest  negative  value  of  the  spectral  fluxes  Il(K)  varies 
by  two  orders  of  magnitude  across  ocean  basins.  As  in 
Arbic  et  al.  (2013),  we  find  that  filtering  model  output 
diminishes  the  deepest  negative  nKE,i(^)  values  and  shifts 
the  model  nKE,i(^)  curves  toward  lower  wavenumbers. 


The  filtered  HYCOM  nKE,i(^)  curves  lie  closer  to  the 
AVISO  curves  than  the  unadulterated  HY COM  curves 
do,  especially  when  both  spatial  and  temporal  filtering 
are  employed.  Exact  matches  between  filtered  HYCOM 
and  AVISO  curves  cannot  be  expected,  not  only  be¬ 
cause  of  model  limitations  but  also  because  our  filtering 
procedure  mimics  the  AVISO  construction  procedure 
only  qualitatively. 

The  negative  nKE,i(^)  values  seen  across  the  entire 
frequency  range  in  the  unadulterated  HYCOM  results 
of  Figs,  llb-f  are  consistent  with  QG  behaviors  and 
demonstrate  that  in  those  regions  nonlinearities  drive 
a  kinetic  energy  cascade  toward  lower  frequencies.  As 
noted  in  section  4,  the  removal  of  temporal  trends  and 
means  is  crucial  for  obtaining  this  consistency  with  QG 
results  across  the  entire  plotted  frequency  range.  As  in 
the  QG  results,  significant  Hke,!  values  take  place  over 
about  two  decades  in  frequency,  but  only  about  one 
decade  in  wavenumber.  In  the  midlatitude  southeast 
Pacific  (Fig.  11a),  nKE,i(^)  unadulterated  HYCOM  is 
negative  for  higher  frequencies  but  positive  for  lower 
frequencies,  in  contrast  to  HY  COM  behaviors  in  other 
regions  and  to  QG  behaviors.  An  explanation  for  the 
positive  nKE,i(<^)  values  at  low  frequencies  in  the 
HYCOM  midlatitude  southeast  Pacific  remains  to  be 
found.  We  do  note  that  in  the  same  region  the  AVISO 
nKE,i(<^)  values  are  also  positive  over  a  similar  range  of 
(lower)  frequencies.  In  like  manner  to  the  wavenumber 
domain  behaviors,  filtering  diminishes  the  negative 
nKE,i(<^)  lobes  in  the  HYCOM  results  and  shifts 
nKE,i(^)  curves  toward  lower  frequencies.  The  AVISO 
spectral  fluxes  nKE,i(^)  ^l^o  tend  to  be  diminished  and 
shifted  toward  lower  frequencies  relative  to  HYCOM. 
In  general,  the  filtered  HYCOM  fluxes  lie  closer  to  the 
AVISO  fluxes  than  do  the  unadulterated  HYCOM 
fluxes.  In  contrast  to  the  AVISO  wavenumber  fluxes 
nKE,i(^)’  clearly  identifiable  negative  lobes  are  not  seen 
in  the  AVISO  nKE,i(^)  curves  across  all  six  regions.  The 
discrepancies  between  the  signs  of  AVISO  and  un¬ 
adulterated  HYCOM  frequency  domain  spectral  fluxes 
could  be  due  to  important  physics  missing  from  HYCOM 
or  to  the  filtered  nature  of  the  AVISO  product.  Filtering 
shifts  the  HYCOM  nKE,i(^)  values  from  negative  to 
positive  over  some  frequencies  and  in  some  regions  (see 
Figs,  lid.  He,  and  especially  Fig.  Ilf).  Finally,  we  note 
that  the  position  of  the  Lj  and  LJU  scales,  respectively, 
relative  to  the  shapes  of  the  Hke,!  curves,  are  similar  in 
the  unadulterated  HYCOM  results  (black  curves  in  Figs. 
10  and  11)  and  QG  results  (black  curves  in  Fig.  9). 

c.  Spectral  transfers 

The  nonlinear  spectral  kinetic  energy  transfers 
^KE,i(^,  from  the  upper  layer  of  the  QG  simulation. 
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Fig.  10.  Nonlinear  surface  ocean  spectral  kinetic  energy  flux  nKE,i(^)  in  wavenumber  K  space,  for  six  regions  shown 
in  Fig.  2,  computed  from  unadulterated  HYCOM  output,  spatially  Altered  HYCOM  output,  temporally  Altered 
HYCOM  output,  HYCOM  output  Altered  in  both  space  and  time,  and  AVISO.  See  text  for  details  of  spatial  and 
temporal  Alters.  Dashed  vertical  lines  denote  the  region-dependent  deformation  wavenumbers  1  rad/L^  (see  text). 


and  from  Agulhas  HYCOM  and  AVISO  ocean  surface  positive  transfers.  The  transfers  are  noisier  than  the 

output,  are  shown  in  Fig.  12.  To  display  the  smaller  neg-  fluxes  because  the  fluxes  are  integrals  of  the  transfers, 

ative  transfers  clearly  while  using  a  symmetric  colorbar,  However,  as  with  the  spectra  and  spectral  fluxes,  it  is  still 

we  choose  a  colorbar  range  that  saturates  for  the  larger  clear  that  the  spectral  transfers  are  wider  in  frequency 
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Fig.  11.  As  in  Fig.  10,  but  for  nonlinear  surface  ocean  spectral  kinetic  energy  fluxes  HkeuC^)  frequency  (co) 
space.  Dashed  vertical  lines  denote  the  region-dependent  frequencies  1  radl(LJU),  where  is  the  deformation 
radius  and  U  is  the  mean  shear  (see  text). 


Space  than  in  wavenumber  space.  Substantial  Tke,!  The  QG  transfers  (Fig.  12a)  are  generally  negative  for 
values  cover  about  one  decade  in  wavenumber,  meaning  higher  wavenumbers  and  positive  for  lower  wave- 
that  a  narrow  range  of  wavenumbers  is  engaged  in  numbers.  At  higher  frequencies  the  QG  transfers  are 
nonlinear  transfer.  The  range  is  wider — two  or  more  negative,  while  at  lower  frequencies  the  sign  of  the 
decades  depending  on  record  duration — for  frequency,  transfer  depends  on  wavenumber.  The  diagonal  cut  in 
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Fig.  12.  Nonlinear  spectral  kinetic  energy  transfer  rKE,i(^’ 
in  (a)  upper  layer  of  idealized  two-layer  QG  simulation  and  from 
Agulhas  region  of  (b)  HYCOM  and  (c)  AVISO  ocean  sur¬ 
face  output.  The  units  are  [(wi  -  W2)^/E^/]/([rad/L^]{rad/[Lj/ 
(W1-W2)]})  in  (a)  and  (nWkg“^)/[(rad  day“^)(radkm“^)]  in 
(b)  and  (c).  Dashed  horizontal  and  vertical  lines  are  as  in  Fig.  4.  The 
X  and  y  axes  are  also  as  in  Fig.  4. 


the  right-hand  side  of  the  QG  spectrum  (Fig.  4a)  is  visible 
in  the  QG  transfers  as  well.  The  HYCOM  and  AVISO 
transfers  (Figs.  12b  and  12c,  respectively)  are  somewhat 
similar  to  the  QG  transfers,  but  are  much  noisier. 

The  QG  wavenumber  domain,  upper-layer  kinetic 
energy  transfer  rKE,i(^)  is  shown  as  a  black  curve  in 
Fig.  13a.  Consistent  with  earlier  work  (e.g.,  Salmon 
1978, 1980;  Hua  and  Haidvogel  1986;  Larichev  and  Held 
1995;  LaCasce  1996),  a  negative  lobe  is  seen  in  rKE,i(^) 
near  the  1/L^  wavenumber  (vertical  dashed  line),  im¬ 
plying  that  nonlinearities  act  as  a  sink  of  kinetic  energy 
near  the  scales  associated  with  baroclinic  instability. 
The  positive  lobe  at  lower  wavenumbers  implies  that 
nonlinearities  deposit  the  energy  removed  at  near-L^ 
scales  into  lower  wavenumbers.  Figure  13a  also  displays 
the  other  terms  in  the  spectral  transfer  budget  [Eq.  (12)]. 
The  nonlinear  upper-layer  kinetic  energy  and  APE 
transfers  are  the  largest  terms  in  the  T{K)  budget  and 
roughly  balance  each  other  out  as  discussed  in  Salmon 
(1980).  The  forcing  and  friction  terms  are  important 


(a)  QG  T(K)  budget,  integrated  over  allco 


Fig.  13.  Wavenumber  domain  spectral  transfers  T{K),  integrated 
over  all  w.  (a)  Upper-layer  spectral  kinetic  energy  transfers 
Fke,i(^) — black  curve — in  idealized  two-layer  QG  simulation, 
shown  against  other  terms  in  the  depth-averaged  spectral  transfer 
budget  [Eq.  (12)].  Dashed  vertical  line  represents  the  deformation 
wavenumber  1  rad/L^.  (b)  Ocean  surface  spectral  kinetic  energy 
transfers  Tke,i(^)  for  Agulhas  region  HYCOM  and  AVISO  out¬ 
put.  Dashed  vertical  line  represents  the  deformation  wavenumber 
1  rad/L^  for  the  Agulhas  region  (see  text). 


secondary  contributors  to  the  spectral  transfer  T{K) 
budget.  The  HYCOM  and  AVISO  wavenumber  domain 
spectral  kinetic  energy  transfers  rKE,i(^)  (Fig*  13b)  are 
similar  to  those  in  the  QG  simulation,  in  that  negative 
lobes  are  seen  near  the  scale  and  positive  lobes  are 
seen  at  lower  wavenumbers. 

Figure  14a  shows  the  frequency  domain  spectral 
transfer  budget  [Eq.  (12)]  for  the  QG  model.  The  neg¬ 
ative  lobe  in  the  upper-layer  nonlinear  spectral  kinetic 
energy  transfer  rKE,i(<^) — black  curve — implies  that 
nonlinearities  extract  energy  near  the  LJU  time  scale 
associated  with  baroclinic  instability  (indicated  by  ver¬ 
tical  dashed  line).  The  nonlinear  terms  deposit  this  en¬ 
ergy  into  lower  frequencies,  as  indicated  by  positive 
black  curve  values  at  lower  frequencies.  To  an  even 
greater  extent  than  in  the  wavenumber  domain  transfers 
T{K),  the  largest  frequency  domain  transfers  T{w)  at 
lower  frequencies  are  the  nonlinear  upper-layer  kinetic 
energy  and  APE  terms,  which  to  the  first  order  balance 
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Fig.  14.  Frequency  domain  spectral  transfers  r(a>),  integrated 
over  all  K.  (a)  Upper-layer  spectral  kinetic  energy  transfers 
7"ke,i(<^) — black  curve — in  idealized  two-layer  QG  simulation, 
shown  against  other  terms  in  the  depth- averaged  spectral  transfer 
budget  [Eq.  (12)].  Dashed  vertical  line  represents  the  frequency 
1  rad/(L^/C/),  where  U  is  the  imposed  mean  shear  ui-U2. 
(b)  Logarithmically  smoothed  ocean  surface  spectral  kinetic  en¬ 
ergy  transfers  Tj^^  i{(x))  for  Agulhas  region  HYCOM  and  AVISO 
output.  Dashed  vertical  line  represents  the  frequency  1  rad/(Lj/t/), 
where  Lj  is  the  deformation  radius  and  U  is  the  estimated  mean 
shear  for  the  Agulhas  region  (see  text). 


each  other  out.  The  forcing  and  friction  terms  are  impor¬ 
tant,  but  only  about  half  as  large  as  the  nonlinear  kinetic 
and  potential  energy  transfers.  Thus,  at  least  in  our  simple 
QG  simulation,  nonlinearities  are  the  most  important 
terms  in  the  maintenance  of  low-frequency  energy. 

Figure  14b  shows  Agulhas  HYCOM  and  AVISO 
7kE,i(<^)  curves,  which  have  been  logarithmically 
smoothed  here  for  display  purposes  (no  smoothing  of  T 
has  been  done  elsewhere  in  this  paper).  At  low  fre¬ 
quencies  the  HYCOM  and  AVISO  rKE,i(<^)  curves  are 
predominantly  positive,  as  in  the  QG  results,  albeit 
noisy.  Equation  (20)  demonstrates  that  any  nKE,i(<^) 
curve  having  predominantly  negative  values  at  higher 
frequencies — such  as  occurs  in  the  Agulhas  region 
for  both  HYCOM  and  AVISO — must  be  associated 
with  predominantly  negative  spectral  transfer  values  at 
high  frequencies.  A  negative  lobe  can  be  seen  in  the 
HYCOM  rKE,i(^)  curve  on  either  side  of  the  baroclinic 


instability  time  scale  LJU  (vertical  dashed  line),  as  in 
the  QG  results.  The  AVISO  transfers  are  predominantly 
negative  at  higher  frequencies,  though  noisy,  and  shifted 
to  lower  frequencies  relative  to  HYCOM,  as  in  the 
spectral  flux  results. 


6.  Summary  and  discussion 

Following  our  earlier  examination  of  geostrophic 
turbulence  in  the  frequency  domain  (ASFMRS),  this 
paper  has  examined  geostrophic  turbulence  in  the 
frequency-wavenumber  domain.  We  are  motivated 
by  Penduff  et  al.  (2011)  and  related  studies,  which 
point  to  the  potential  of  oceanic  mesoscale  eddies  to 
drive  significant  low-frequency  oceanic  variability.  The 
frequency-wavenumber  diagnostics  developed  here  and 
in  ASFMRS  represent  a  relatively  little-known  tool  for 
studying  low-frequency  oceanic  variability  [to  our 
knowledge,  the  only  previous  usage  of  such  diagnostics 
is  in  the  atmospheric  literature,  e.g.,  Sheng  and  Hayashi 
(1990a,b)].  Frequency-wavenumber  spectra,  spectral 
fluxes,  and  spectral  transfers  from  an  idealized  forced- 
dissipated  QG  turbulence  simulation,  a  realistic  eddying 
HYCOM  simulation,  and  AVISO  gridded  satellite  al¬ 
timeter  maps  exhibit  several  similarities.  The  range  of 
frequencies  that  are  energetic  and  participate  strongly  in 
nonlinear  interactions  is  greater  than  the  range  of  en¬ 
ergetic  and  active  wavenumbers.  Consideration  of  joint 
frequency-wavenumber  spectral  fluxes  yields  the  im¬ 
portant  constraint  that  the  frequency  domain  spectral 
flux  at  the  lowest  frequency  must  be  zero  for  geostrophic 
flows.  Application  of  this  constraint  gives  HYCOM 
frequency  domain  fluxes  and  transfers  that  are  more 
robust,  and  more  consistent  with  QG  model  behaviors, 
than  the  frequency  domain  spectral  fluxes  computed 
from  realistic  models  in  ASFMRS.  In  five  of  the  six  re¬ 
gions  examined  in  this  paper,  the  realistic  HYCOM 
simulation  exhibits  a  kinetic  energy  cascade  toward 
lower  frequencies  alongside  the  well-known  inverse 
cascade  toward  lower  wavenumbers,  as  in  the  QG 
model.  AVISO  spectral  fluxes  show  a  cascade  to  lower 
frequencies  in  some  regions,  but  not  in  others.  However, 
the  discrepancies  of  the  AVISO  results  with  QG  and 
HYCOM  results  are  likely  due  at  least  in  part  to  the 
filtering  inherent  in  the  creation  of  the  AVISO  product. 
Exercises  performed  with  HY COM  output  indicate  that 
such  filtering  distorts  the  computed  frequency  domain 
spectral  fluxes  substantially.  In  the  QG  simulations, 
spectral  transfer  budgets  demonstrate  that  nonlinear 
transfers  of  kinetic  and  potential  energy  are  the  primary 
contributors  to  the  maintenance  of  low-frequency,  low- 
wavenumber  energy,  with  forcing  and  friction  playing 
important  but  secondary  roles.  Thus,  in  the  QG  model  at 
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least,  low-frequency  variability  is  more  due  to  intrinsic 
nonlinearities  than  to  forcing. 

Low-frequency  energy  is  associated  with  low  wave- 
numbers,  and  vice  versa,  though  not  in  a  simple  way. 
Spectral  slopes  in  wavenumber  and  frequency  space  are 
not  equal,  as  they  would  be  if  Taylor  (1938)’s  hypothesis, 
which  relates  frequency  to  wavenumber  spectra,  held 
precisely.  The  limitation  of  the  Taylor  (1938)  hypothesis  in 
our  results  is  consistent  with  discussions  in  Hayashi  and 
Golder  (1977),  who  found  no  simple  relationship  between 
the  frequency  and  wavenumber  spectra  of  kinetic  energy 
in  realistic  atmospheric  models  and  atmospheric  datasets. 
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